CLAIMS 

1.-70. (CANCELLED) 
71. (NEW) A method of estimating blood perfusion indices, for a region of interest (ROI) 
by operating a computer program on intensity data input to a computer comprising: 

a. providing a contrast agent to an artery, wherein said artery is upstream from a 
region of interest (ROI), wherein said ROI comprises tissue; 

b. measuring intensity data to provide a tissue contrast agent curve C(i) in said 
tissue, wherein said C(t) is a known intensity profile measured against time; 

c. measuring an arterial input function AFT a (t) in a vessel leading to said ROI, 
wherein said AFT a (t) is known intensity profile measured against time; 

d. determining a tissue blood flow F t and a tissue impulse residue function R(t) by 
deconvolving C(t) = (F t /k H )AIF a (t)®R(t), wherein k H =(l-H a )l(l-H t ) is a 
correction constant using different values of an arterial hemotaocrit H a and a 
tissue hematocrit H t , wherein said H a has a value of about 0.45, wherein said H t 
has a value of about 0.25, wherein when said AFT a (t) in a major artery comprises 
no delay and no dispersion, said C(t) can be derived by a convolution of said 
AIF a (t) and said R(t); 

e. determine a blood flow rate F t and a tissue impulse residue function R^t) using a 
said C(t) and said AIF t by deconvolving the relation C(t)= (F t /k H )AIF t (t) ® Re(t), 
wherein ku is a hermocrit correction constant having a known value; 

f. measuring an arterial input function AIF a (t) in said artery at a position that is 
upstream from said tissue; 

g. simulating a tissue contrast agent concentration curve C s (t) using a vasucular 
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transport function h a {t), and a tissue transport function h s {t), wherein said h s {t) is 
given by a Gaussian probability density function, wherein said 

h s (t)=— (t-t 2 ) a2 e-°- t2Va2 for t s> t 2 , and said h s (t)= 0 for t < t 2 , and 
Ai 

h s (t)=— (t -t 2 ) a2 e~ (t ~ t2), ° 2 for t ^ 0, wherein said a 2 and a 2 are parameters 

5 related to mean transit time and dispersion of said h s (t) 9 wherein (a 2 * a 2 ) is a 

peak rise time, wherein a 2 (l + a 2 ) is a mean transit time, wherein for said a 2 , 
wherein when said ct 2 = 0, said h s (i) = — e ~ {t ~ h)lal for t> t 2 , and said h s (t) = 0 for 

t < t 2 , wherein l/a 2 is a peak height, wherein t 2 + a 2 is a mean transist time (r) in 
said tissue; 

10 h. using an adiabatic approximation to a homogeneity model in said tissue to 

provide a tissue input residual function R (t, r), wherein for t =s r, said R (t, r) = 
1, wherein fro t > r, said R (t, r) = Ee~ k(t ~ r \ wherein E is an extraction fraction 
of said contrast agent in said blood that leaks out of a vessel into said tissue, 
wherein k is a clearance rate at which said contrast agent diffuses back into said 

15 blood and leaves said tissue, wherein said k = E*F t /V e , wherein said V e us a 

volume fraction of extravascular space and extracellular space in said tissue, 
wherein said E comprises a value between 0 and 1, wherein said V e comprises a 
value between 0 and 1 ; 
i. providing in average tissue input residue function ^(0, wherein said R s (t) 

20 comprises a Gaussian probability density function h s (f) of a contrast agent, 
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t t 

wherein R s (t) = 1 - f h s (r)dr+ Ee~ kt j h s (r)e kT dr ; 

0 0 

j. determining a simulated contrast agent curve C s (t), wherein said C s (t) is 
reperesented by C s (t) = (F t /k H )AIF t (t)®R s (t) = (F, lk H )] AIF t {x)R s (t - x)dx\ 

0 

k. using an iterative least squares method to fit the said simulated C s (t) to said 
5 measured tissue curve C(t) with four adjustable parameters; 

1. using a model-free deconvolution method to estimate starting values of said 
adjustable parameters for said iterative least squares fitting process in order to 
derive optimized values of said adjustable parameters; and 
m. calculating perfusion indices from the said optimized values of the adjustable 
10 parameters. 

72. (NEW) The method of claim 70, wherein said intensity data is generated by 
administering a contrast agent to a body lumen of a body during a dynamic 
imaging scan, wherein said body lumen comprises an artery or a vein, wherein 

15 an image response from said contrast agent is recorded to computer data storage 

in a computer. 

73. (NEW) The method of claim 70, wherein said C(t) is a temporal concentration of 
said contrast agent obtained from said intensity data, wherein said intensity data 

20 comprises contrast images sequentially acquired from a region in a body, 

whereby said contrast agent concentration is plotted versus time. 
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74. (NEW) The method of claim 70,wherein said AIF a is scaled upward according to 
a venous input function (VIF), wherein said VIF is based on a measured contrast 
intensity profile from a vein draining from said ROI. 
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